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Rainfall retrievals from polarimetric radar measurements require the knowledge of 
four fundamental rain microstructure parameters, namely, drop size distribution, 
drop shape distribution, canting angles and drop fall velocities. Some recent 
measurements of all four parameters in natural rain are summarized in [1]. In this 
paper, we perform an in-depth analysis of two events, using two co-located 2D 
video disdrometers (2DVD; see [2]) both with high calibration accuracy, and a C- 
band polarimetric radar [3], located 15 km away. The two events, which occurred 7 
days apart (on the 18th and the 25th of Dec 2009), had moderate-to-intense 
rainfall rates, but the second event had an embedded convection line within the 
storm. The line had passed over the 2DVD site, thus enabling the shapes and fall 
velocities to be determined as the line crossed the site. The first event was also 
captured in a similar manner by both the 2DVDs as well as the C-band radar. 

Drop fall velocity measurements for, say, the 3 mm drops show noticeable 
differences between the two events. Whereas for the first event, the velocity 
distribution showed a narrow and symmetric distribution, with a mode at the 
expected value (7.95 m/s, as given by the formula in [4]), the second event 
produced a wider distribution with a significant skewness towards lower velocities 
(although its mode too was close to the expected value). Moreover, the ‘slower’ 3 
mm drops in the second event occurred when the convection line was directly over 
the 2DVD site (03:35-03:45 utc), and not before nor after. A similar trend was 
observed in terms of the horizontal dimensions of the 3 mm drops, i.e. large 
fluctuations during the same time period, but not outside the period. Vertical 
dimensions of the drops also fluctuated but not to the same extent. Interestingly, 
the horizontal dimensions tended towards larger values during the 10-minute 
period, implying an increase in drop oblateness, which in turn indicates the 
possibility of the ‘horizontal’ mode oscillation, one of the three fundamental modes 
of drop oscillations [5], albeit the most difficult one to excite. 

In order to dismiss the possibility of surface layer effects, the C-band polarimetric 
radar data were processed for this event, and analyzed in terms of the ‘effective 
beta’ method [6]. This method, although susceptible to radar calibration errors, can 
be used qualitatively to identify areas which do not conform to the polarimetric 
radar algorithms built on standard bulk assumptions on drop shapes. Application of 
this method to the PPI scans taken over a one-hour time period related to the 
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second event has shown higher than normal beta-values in some areas of the 
convection line, the high values corroborating the notion of horizontal mode 
oscillations occurring in those regions. 

We will present the fall velocity measurements and drop shapes for both events 
from the two 2DVDs, as well as results of the corresponding C-band data 
analyses. 
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1. INTRODUCTION 

The estimation of rainfall rates from polarimetric 
radar measurements is implicitly based on certain 
assumptions regarding drop shapes, orientations 
and drop fall velocities within the radar pulse 
volume. Typically it is assumed that: (i) the drop 
axis ratio versus drop equivalent diameter (D eq ) 
follows a single, monotonic variation (eg. Brandes 
et al., 2002), (ii) the drop canting angles have a 
Gaussian distribution with a mean of zero 
degrees and a standard deviation of around 7-8 
degrees, and (iii) the fall velocity has a unique 
variation with D eq , for example, the formula given 
in Atlas et al. (1973) at sea level. Polarimetric 
radar measurements are then used to retrieve the 
drop size distribution (DSD) within the pulse 
volume assuming such ‘bulk’ assumptions, and 
the rainfall rate within the pulse volume is 
estimated. 

While the above approach results in reasonably 
accurate estimates or rainrate (see, for example, 
Ryzhkov et al. 2005, Bringi et al. 2011), there 
have also been a few studies which show that 
such bulk assumptions may not apply universally. 
For example, Gorgucci et al. (2000) use Z h , Z dr , 
and K dp to estimate the ‘effective’ drop shapes 
and subsequently modify the rainrate retrieval 
algorithms using the so-called ‘effective beta- 
method’. Another study (Thurai et al. 2009) 
demonstrated that scattering calculations based 
on shape (and size), and orientation of individual 
drops determined from a well-calibrated 2d-video 
disdrometer (2DVD) gave closer agreement with 
polarimetric radar measurements (Z h , Z dr , K dp as 
well as p hv ) than those based on the measured 
DSDs and using bulk assumptions. Indeed the 
study high-lighted the possibility of ‘significant 
deviation’ from the mean shape assumptions in 
one event study. 

It is well established that the 2DVD (Shonhuber et 
al., 2008) if calibrated accurately can measure 
shape, orientation and fall velocity of individual 


drops transiting its sensor area. Such 
measurements have been reported previously in 
Thurai and Bringi, (2008), and Huang et al. 
(2008). In this paper, we report 2 case-events 
which occurred within 7 days of each-other, one 
conforming to the ‘bulk-assumptions’ regarding 
drop axis ratios and fall velocities and the other 
showing significant deviations within one part of 
the storm. Measurements were made using two 
collocated (i.e. a few meters apart), and 
accurately calibrated 2DVDs as well as a C-band 
polarimatric radar (ARMOR, see Petersen et al. 
2007) located 15 km from the 2DVD site. The 
measurements were part of an on-going long- 
term campaign in Huntsville, Alabama. 


2. THE TWO EVENTS 

The two events which occurred on 18 December 
2009 and 25 December 2009 were of long 
duration lasting several hours. Fig. 1(a) and 1(b) 
show composite radar images of the two events 
around the Huntsville area. The white cross in 
both panels marks the 2DVD location. In both 
cases, the event had high reflectivities, but for 
case 2, a well-defined embedded thin line of 
convection can be seen crossing the Huntsville 
area. Both events had relatively high rainfall 
accumulations. Table 1 compares the daily totals 
from the two 2DVD data and from a Geonor rain- 
gage, also collocated. 


Table 1 : Total rainfall (mm) from the two 2DVDs 
and from the collocated Geonor 


Date 

SN16 

SN25 

Geonor 

Event 1 : 

18 Dec 2009: 

35.5 

>33. 2* 

34.5 

Event 2: 

25 Dec 2009: 

26.2 

24.4 

25.4 


( T some loss of data for several minutes) 





Fig 1 : Composite radar images of the two events 
considered in this study. The white cross shows 
the location of the two 2DVDs. The ARMOR radar 
is situated 15 km away in the south-west 
direction. Note event 2 has a line convection 
crossing the 2DVD location. 


3. 2DVD MEASUREMENTS 

As mentioned above, the 2DVD measures the 
drop fall speed (i.e. the vertical velocity) for each 
drop falling within the sensor area, as well as the 
shape, size and orientation. Fig. 2(a) and 2(b) 
shows the distribution of the measured vertical 
velocity for all the 3 mm drops (to be precise 3 ± 
0.1 mm) for both events. In both panels, the blue 
curve represents the distribution determined from 
one of the instruments (the low-profile 2DVD, also 


labeled SN-16) and the red curve represents the 
distribution determined from the other instrument 
(the compact 2DVD, labeled SN-25). In both 
cases the red and the blue curves show good 
agreement (thus providing confidence in the 
2DVD measurements) but the two events show 
different characteristics in terms of their 
probability distributions. Whereas for event 1, the 
distributions are narrow, symmetric and have a 
mode close to the expected fall velocity for D eq =3 
mm (7.95 m/s), the second event shows a wider 
distribution with a noticeable negative skewness. 
While both events likely had different wind 
conditions with different up-down drafts, and 
turbulence etc., the skewness in Fig. 2(b) cannot 
be explained easily. Moreover, as will be shown 
later, the second event also shows significant 
deviation from our standard ‘mean-shape’ during 
a certain period of the event. 



Fig. 2: Distribution of measured velocities of all 
the 3 mm drops, from the two 2DVD 
measurements for event 1 (top) and event 2 
(bottom). 



We now investigate the fall velocity of each of the 
individual 3 mm drops captured by the two 
2DVDs for event 2. The upper panels in Fig. 3 
show these velocities plotted as time series as 
the event passed over the 2DVD site. In both sets 
of measurements, the same observations can be 
made, that is, larger velocity fluctuations between 
03:30 and 03:50 UTC, with a significant 
proportion of the drops having lower than the 
expected velocities. Outside this time range, the 3 
mm drop velocities show much less variation and 
are centered around the 7.95 m/s expected value. 
Clearly, the negative skewness in Fig. 2 arises 
primarily due to the ‘slower’ 3 mm drops captured 
between 03:30 and 03:50 UTC. This is significant, 
and cannot be dismissed as being due to any 
instrument ‘calibration problems’. Note the 
velocities represent the actual fall speed, with no 
contribution from the wind-induced horizontal 
component. 


Using the unique image processing capability of 
the 2DVD, i.e. the contoured images of each 
drop, both datasets were used to determine the 
maximum horizontal and vertical dimensions of 
each drop, after removing the effects due to drop 
orientations as well as the drop horizontal velocity 
components (Fig. 3, lower panels). Again, both 
sets of measurements show similar 
features/characteristics, i.e. between 03:30 and 
03:50 UTC, there appear large variations in the 
drop dimensions as compared with other times. 
The larger variation is particularly noticeable for 
the horizontal dimension, with a significant 
number of drops showing larger than expected 
values (the mean horizontal dimension is around 
3.2 to 3.3 mm). These larger horizontal 
dimensions will result in lower axis ratios, i.e. an 
increase in oblateness , a somewhat surprising 
and unusual finding. In the next section, we 
examine the C-band polarimetric radar 
observations. 
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Fig. 3: The measured fall velocity of each individual 3 mm drops, as time series (top panels) and 
their corresponding horizontal and vertical dimensions (red and green respectively) measured by 
the SN16 (left panels) and SN-25 (right panels) for event 2. Fluctuations are noticeable between 

03:30 and 03:45 UTC. See main text. 




4. RADAR OBSERVATIONS (EVENT 2) 

The composite image in Fig. 1b for event 2 
shows that the convective line crossed the 
2DVD site, and the time sequence of Fig. 3 
demonstrates that it crossed the Huntsville 
area between 03:00 and 04:00 UTC. The 
ARMOR radar made routine observations 
during the day, and in Fig. 4 we show a set of 
panels of 1.3 deg elevation PPI scans taken 
at 03:05, 03:40 and 03:55 UTC. The upper 
panels show the reflectivity Z h , the middle 
panels show Z dr and the lower panels show 
K dp . The organized convective line can be 
seen to move from south-west to north-east, 
and at 03:40, the line can be seen to lie 
directly above the 2DVD site. It is around this 
time that the 2DVD measurements (Fig. 3) 
show the larger velocity fluctuations and the 
large axis ratio variations, resulting in ‘slower 
drops’ with higher oblateness. At 03:55 UTC, 
the line is to the north-east of the 2DVD site 
and is seen to be fragmenting and/or 
decaying. In the 03:05 case, i.e. when the line 
was well-defined and well-organized, the Z dr 
correction procedure failed to correct for the 
total differential attenuation beyond the 
convective line. Since the correction 
procedures are based on our mean drop 
shapes (and big drop regime assumptions), 
the difficulty in correcting for differential 
attenuation indicates an increase in drop 
oblateness in the region of the line 
convection. The implications for proper 
attenuation correction such instances are 
interesting to consider from the perspective of 
operational mapping of rainfall characteristics 
using C-band polarimetric radar. 


5. INTERPRETATION OF OBSERVATIONS 

Ancillary meteorological data for Huntsville 
indicates that this event was associated with 
relatively high winds and gusty conditions 
near the time of convective line passage, 
thence there is a possibility that surface layer 
(i.e. close to ground level) turbulence could 
have facilitated unusual drop oscillation 
modes, such as mixed and/or horizontal 
oscillation modes (Beard et al. 2010) which 


result in increased drop axis ratios, and a 
possible corresponding reduction in fall 
velocity. In the literature, the effect of drop 
axis ratio variations on the fall velocity has not 
been addressed in a quantitative manner due 
to the complexity, but some interesting 
discussion can be found in section 5(a) of 
Beard (1976) regarding the interpretation of 
some previously published results, for 
example by Jones (1959), which implies that 
turbulence may cause drop shape distortion. 
Interestingly, Jones’ data showed the axis 
ratios of the 4 mm drops ranging from 0.5 to 
1.3 determined from 2000 photographs taken 
at ground level during turbulent events, 
compared with our results for the 3 mm drops 
ranging from 0.6 to 1.2 during the line 
convection period. 

The question remains whether the mixed 
mode oscillations (including the horizontal 
mode) associated with the line convection in 
event 2 occurs only within the surface layer 
near the 2DVD or at higher heights too. To 
address this, we analyze the radar data in 
Fig. 4. The higher than normal axis ratios 
within the line convection should be reflected 
in higher Z dr for a given DSD. On the other 
hand, the retrieval of the DSD, and in 
particular, the mean volume diameter, D 0 , 
from Z dr is based on an assumed axis ratio 
dependence on D eq . To overcome this 
‘circular ambiguity’, we have tried to apply the 
so-called ‘effective-beta method’ to the PPI 
scans. This method involves the calculation of 
the 3-effective parameter (3eff) which 
represents the approximate slope of the curve 
representing the axis-ratio versus D eq 
variation. The parameter 3eff is calculated 
using Z h , Z dr and K dp , and is mostly 
independent of the DSD within the pulse 
volume. The calculated values can enable 
one to assess whether the radar 
observations conform to the ‘bulk 
assumptions’ regarding drop axis ratios. 
However, for the method to be applicable, 
both Z h and Z dr need to be accurately 
calibrated and properly corrected for 
attenuation, and, moreover, the expected 
values of (3eff) depend on the exact equation 
used to derive it. 
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Fig. 4: PPI scans of Zh (1 st row), Zdr (2 nd row) and Kdp (3 rd row), taken at 03:05, 03:40 and 03:50 
UTC. The 2DVD site is marked with a star (azimuth 52°, range 15 km). Note Z dr suffers severe 
differential attenuation beyond the convection line, particularly for the 03:05 case, even after 
applying attenuation correction procedures based on our mean drop shapes, implying the 
presence of more oblate drops than normal. Small amounts of cross-coupling effects may also be 

possible, perhaps. 


The calibration of Z h and Z dr were established 
by comparing the radar data extracted over 
the disdrometer site with corresponding 2DVD 
measurements collected well after the line of 
convection passed the 2DVD site between 
04:00 and 05:00 UTC. The calibration factors 
were found to be relatively steady throughout 


the hour for both Z h and Z dr . The PPI scans 
given in Fig. 4 represent radar data after 
applying the calibrations. 

The range profiles extracted from the 03:40 
PPI scan along azimuths between 50 to 55 
degrees are shown in Figure 5. At the 2DVD 
site, which is located at 15 km from the radar, 


reflectivity values of 50-55 dBZ can be 
observed, together with Z dr of around 4 dB and 
K dp of 4 to 5 deg/km. In this region, the 
calculated values of p e ff (using an equation 
based on DSD model based simulations) are 
also significantly higher, as compared with say 
those at 17-25 km range. The higher values of 
peff seen at the 2DVD site indicate a steeper 
slope of the axis ratio versus D eq variation, 
which is consistent with the presence of more 
oblate drops in the 2DVD measurements (Fig. 
3). However, some caution with respect to the 
calculated (3 e ff values is warranted for regions 
beyond the line of convection, because of the 
lowered Z dr even after correcting for 
differential attenuation, as was mentioned 
earlier. 




Range, km 

Fig. 5: Range profiles taken along 50-55 deg 
azimuth of Z h (top), Z dr (2 nd ), K dp (3 rd ) and the 
calculated (3 e ff (lowest), after attenuation 
correction and calibration. The 2DVDs are 
located at 15 km range. 


The peff values calculated for the PPI scan at 
03:40 are shown in Fig. 6. Some smoothing 
has been applied in order to better identify 
regions of ‘abnormal’ (3 e ff values. Additionally, 
the (B e ff values were only calculated in regions 
with Z h > 35 dBZ, Z dr > 0.3 dB and K dp > 0.3 
deg/km. Most regions within the line of 
convection met these criteria, as did a number 
of other isolated areas. Within the convective 
line, small areas of relatively high p e ff values (> 
0.7) can be observed very close to the 2DVD 
site (marked with a black star). Further to the 
south, a larger area of higher p e ff values can 
also be seen whereas in other regions within 
the line convection, more moderate p e tf values 
or 0.05 to 0.06 are obtained. 
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Fig. 6: The calculated p e ff for the scan at 03:40 
UTC. Red regions are a qualitative indication 
of higher slope for axis ratio versus D eq 
variation. 


6. FURTHER COMMENTS 

As mentioned above, the p e tf values in Fig. 
5(d) and 6 depend on the specific equation 
used, but in relative terms, the higher values 
(i.e. the regions in red color) do indicate that 
the slope of the axis ratios versus D eq is larger 
than the ‘mean’ variation. The presence of 
more oblate drops (perhaps with mixed-mode 
oscillations) is likely in these regions, and the 
2DVD measurements in Fig. 3 support this, 
albeit at ground level. 


It should be stressed here that this particular 
event (i.e. event 2 on 25 Dec 2009) is a 
somewhat unusual case, which was 
investigated because of the negative 
skewness in the fall velocity distributions (Fig. 
2) from the 2DVD measurements. In several 
other cases which have been examined to 
date, no evidence was found to indicate 
mixed-mode drop oscillations resulting in 
higher axis ratios and lower velocities. In the 
literature, evidence of unusual drop 
oscillations (especially relating to horizontal 
mode) is relatively scarce, but not non- 
existent. In one particular study, some 
inferences regarding horizontal oscillations for 
large drops were made (Jameson and Durden, 
1996) from nadir-pointing radar measurements 
in a convective tropical storm. Relatively high 
LDR measurements were seen in certain 
sections of the convective storm which were 
attributed to the possibility of ‘horizontal 
distortions’ particularly for large drops. Further, 
the authors state that ‘collisions between large 
and small drops are adequate to produce and 
sustain the horizontal oscillations...’. It is 
conceivable that such a scenario also applies 
to our case, i.e. to some regions within the line 
of convection in Fig. 5, 6, part of which was 
fortuitously captured by the two 2DVDs. 
Another case (08 April 2010, not reported 
here) has also indicated similar characteristics 
during our preliminary analysis. The fact that 
there are two collocated 2DVDs greatly 
enhances our confidence in these 
measurements (provided there is agreement 
between the two). Such data, together with 
simultaneous ARMOR observations will form 
part of an on-going study to identify cases 
where significant deviations from our mean 
shapes and the expected fall velocities seem 
to occur. 


7. SUMMARY 

Two events which occurred 7 days apart in 
Huntsville, Alabama, have been investigated 
using two collocated 2DVDs, both at ground 
level, as well as simultaneous observations 
from the ARMOR C-band polarimetric radar. 
For each event drop fall velocities and shapes 
were examined from the 2DVD 
measurements, with specific focus placed on 
3 mm diameter drops. The first event - on 18 
Dec 2009 - showed the expected fall 
velocities; for example, the 3 mm drops had 


velocities which were symmetrically 
distributed, with a mode at around 7.9 to 8 
m/s. The distribution was narrow. The second 
event - on 25 Dec 2009 - showed a skewed 
distribution, with a negative skewness with 
significant number of drops having lower fall 
velocities. The event was associated with an 
embedded line of convection, high winds and 
gusty conditions. Time series of the 3 mm drop 
fall velocity measurements showed that these 
‘slow’ drops were detected only during 
passage of this line. The drop maximum 
dimensions were also examined (after 
removing drop canting effects as well as drop 
horizontal velocities). Large fluctuations in the 
horizontal and vertical dimensions were 
observed during the same time period; though 
fluctuations in the horizontal dimensions were 
larger, and were observed to trend towards 
higher values, indicating the presence of 
horizontal mode oscillations, or at least mixed 
mode oscillations, which include horizontal 
mode oscillations. 

Simultaneous radar observations from the 
ARMOR radar were also analyzed for altitudes 
above the surface. Here the so-called (3- 
effective method was applied, which enables 
one to qualitatively identify areas which do not 
conform to the standard bulk assumptions 
often used in polarimetric radar algorithms 
designed to diagnose rain characteristics. The 
results indicated higher than normal p eff in 
some regions within the convective line, 
including a small region directly above the 
2DVD site. The implication is that horizontal 
mode drop oscillations (or at least mixed-mode 
oscillations) were occurring in some regions 
within the convective line, not only at ground 
level but also aloft. The difficulty in correcting 
for differential attenuation, especially when the 
convection line was highly organized (at 
03:05), also corroborates the notion of drop 
horizontal oscillations occurring within this 
region. Studies are on-going to identify other 
cases which show significant deviations from 
our mean shapes and the expected fall 
velocities, observed from both 2DVDs, 
accompanied by the C-band polarimetric radar 
observations. 
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